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ABSTRACT 

The  formulation  of  a  one-dimensional  analytical  model 
representing  steady  state  operation  of  a  vertically  oriented 
gas  loaded  variable  conductance  heat  pipe  is  presented.   The 
model  includes  not  only  the  effects  of  binary  mass  diffusion 
and  axial  pipe-wall  conduction,  but  also  gravitational 
effects  which  exist  at  the  vapor-gas  interface  region  when 
the  molecular  weights  of  the  working  fluid  and  the  noncon- 
densible  gas  differ  significantly.   Analytical  expressions 
for  the  transport  processes  (mass,  momentum,  and  energy) 
along  with  equations  of  state  result  in  a  system  of  ten 
equations  which  describe  steady  state  operation  of  the  device 
Combining  these  expressions  yields  a  system  of  two  simul- 
taneous highly  non-linear  differential  equations.   These 
equations  are  non-dimensionalized  and  set  up  for  numerical 
integration  using  a  variable  step  fourth  order  Runge  Kutta 
method  as  part  of  IBM's  Continuous  System  Modeling  Program 
(CSMP) . 
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I.       INTRODUCTION 

A.  BACKGROUND 

The  heat  pipe  is  a  heat  transfer  device  capable  of  trans- 
ferring large  quantities  of  heat  under  nearly  isothermal 
operating  conditions.   This  device  is  basically  a  liquid-vapor 
two  phase  system  which  transfers  heat  by  evaporation  and 
condensation  of  a  working  fluid  at  opposing  ends  of  a  tube. 
There  are  no  moving  parts  in  a  heat  pipe.   To  sustain  opera- 
tion the  condensed  working  fluid  is  transferred  from  the 
condenser  to  the  evaporator  section  of  the  device  by  capillary 
pumping  action  through  a  wick  structure  lining  the  heat  pipe 
inner  wall. 

B.  CONVENTIONAL  HEAT  PIPE  PRINCIPLES 

The  conventional  heat  pipe  in  its  simplest  form  consists 
of  a  hollow  tube  with  wicking  material  lining  its  inner  sur- 
face as  shown  in  Figure  1.   This  device  is  evacuated,  then 
filled  with  an  appropriate  amount  of  a  working  fluid. 

In  operation,  heat  at  a  given  rate  from  a  source  is  ab- 
sorbed at  one  end  of  the  pipe  evaporating  the  fluid  from  the 
wick.   The  vapor  then  flows  along  the  pipe  to  the  condenser 
end  where  heat  is  rejected  to  a  sink  and  the  vapor  is  condensed 
into  the  wick.   Circulation  of  the  working  fluid  is  maintained 
by  a  capillary  head  in  the  wick  which  must  exactly  balance 
pressure  drops  in  the  liquid  and  vapor,  as  well  as  any  body 
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forces  present.  The  heat  load  alone  controls  the  satura- 
tion temperature  and  pressure  of  the  working  fluid,  which 
become  the  operating  temperature  and  pressure  of  the  pipe. 

Since  the  transfer  of  heat  is  effected  by  adding  and 
removing  the  heat  of  vaporization  at  saturation  conditions, 
the  heat  pipe  is  nearly  isothermal  along  its  length.   How- 
ever, a  change  in  the  heat  load  will  automatically  change 
the  operating  temperature  (and  pressure)  of  the  pipe,  given 
a  constant  sink  temperature.   This  characteristic  is  unde- 
sirable for  many  applications  where  a  constant  source  tem- 
perature is  required. 

C.   VARIABLE  CONDUCTANCE  HEAT  PIPE  PRINCIPLES 
1.   General 

If  a  noncondensible  gas  is  introduced  into  a  heat 
pipe  along  with  the  working  fluid,  the  heat  pipe  becomes  a 
passively  controlled  device  which  can  operate  at  various  heat 
loads  while  maintaining  a  more  nearly  constant  source 
temperature. 

During  operation  the  working  fluid  behaves  as  described 
above;  however,  the  noncondensible  gas  is  entrained  in  the 
vapor  flow  and  forms  a  gas  plug  at  the  condenser  end  of  the 
pipe.   This  plug  acts  as  a  diffusion  barrier  to  the  flowing 
vapor  so  that  the  effective  condenser  length,  and  therefore 
condenser  heat  transfer  area,  is  reduced.   For  a  given  heat 
load  and  noncondensible  gas  loading,  the  pipe  will  operate 
at  a  specific  temperature.   If  the  heat  load  is  increased  the 
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saturation  temperature  and  partial  pressure  of  the  working 
fluid,  as  well  as  total  pressure  in  the  pipe,  tend  to  in- 
crease.  In  response,  the  gas  plug  is  compressed  which 
increases  the  effective  condenser  heat  transfer  area,  thus 
reducing  the  system  temperature  rise.   A  thorough  presentation 
on  the  theory  and  design  of  variable  conductance  heat  pipes 
has  been  provided  by  Marcus  [1]. 

2.  Flat  Front  Theory 

Bienert  [2]  developed  the  earliest  model  used  to 
examine  the  operating  characteristics  of  a  gas  loaded  heat 
pipe.   This  flat  front  theory  assumed  that  the  noncondensible 
gas  was  in  the  form  of  a  plug  with  a  sharply  defined  vapor- 
gas  interface  perpendicular  to  the  axis  of  the  pipe.   This 
simple  model  neglected  axial  conduction  in  the  heat  pipe  wall 
and  the  diffusion  of  vapor  into  the  stagnant  gas  plug.   This 
early  model  provided  much  insight  into  the  operating  charac- 
teristics of  gas  loaded  heat  pipes.   However,  experimental 
results  [1]  have* shown  the  model  to  be  inaccurate,  primarily 
because  it  neglects  axial  conduction  and  diffusion. 

3.  Diffuse  Front  Theory 

A  more  accurate  model  of  the  gas  loaded  heat  pipe  was 
provided  by  Edwards  and  Marcus  [3].   The  model,  referred  to 
as  the  'diffuse  front  theory',  includes  effects  of  binary 
diffusion  between  the  working  fluid  and  the  noncondensible  gas 
as  well  as  axial  conduction  in  the  pipe  wall.   Marcus  reported 
[1]  good  correlation  between  experimental  results  and  a  com- 
puter program  based  upon  the  diffuse  front  theory. 
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4.   Gravitational  Effects 

Neither  the  flat  front  theory  nor  the  diffuse  .front 
theory  included  body  force  effects  on  the  vapor-gas  inter- 
face region.   This  simplification  was  probably  made  because 
much  of  the  early  work  with  heat  pipes  was  related  to  space 
program  applications. 

Kelleher  and  Batts  [4]  reported  convincing  experimen- 
tal results  indicating  that  gravity  effects  can  significantly 
distort  the  axial  temperature  profile  when  the  working  fluid 
and  noncondensible  gas  are  of  significantly  different  molecular 
weights . 

D.   OBJECTIVE 

The  objective  of  this  research  was  to  develop  a  one-dimen- 
sional analytical  model  describing  steady  state  operation  of 
a  gas  loaded  variable  conductance  heat  pipe.   The  model  was 
to  include  not  only  the  effects  of  axial  conduction  and  binary 
mass  diffusion,  but  also  gravitational  effects  on  the  vapor- 
gas  interface  region.   An  additional  purpose  of  this  effort 
was  to  transform  the  analytical  equations  into  a  set  of  non- 
dimensionalized  simultaneous  differential  equations  which 
could  readily  be  solved  numerically  using  IBM's  Continuous 
System  Modeling  Program  III  [5]. 
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II.   FORMULATION  OF  ANALYTIC  EQUATIONS 

In  formulating  the  analytic  equations  describing  steady 
state  heat  pipe  operation  the  differential  control  volume 
in  Figure  2  was  developed  to  represent  a  section  of  the  heat 
pipe  condenser.   As  implied  in  Figure  2,    only  vertically 


X 


X  +  dx 


^  Pi^v,  +  ^T— (p,  V,  )  dx 
^b  b   dx  b  b 


Pa^a^dS^^PaV^^ 


Figure  2.   Differential  Volume  Element 
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oriented  heat  pipes  are  considered  so  that  the  transport 
processes  can  be  treated  as  one  dimensional. 

A.   SIMPLIFYING  ASSUMPTIONS 

The  following  assumptions  were  made  in  formulating  an 
analytical  model  using  the  various  transport  processes: 

1.  The  heat,  mass,  and  momentum  transport  processes  were 
assumed  to  be  one-dimensional  in  the  axial  direction  within 
the  vapor  space.   This  assumption  restricted  this  analytical 
model  to  the  specific  case  of  a  vertically  oriented  heat  pipe. 
In  any  other  orientation,  gravitational  effects  should  be 
expected  to  cause  highly  three  dimensional  spatial  variations 
when  the  molecular  weights  of  the  working  fluid  and  noncon- 
densible  gas  differ  substantially.   This  effect  was  demon- 
strated experimentally  by  Kelleher  and  Batts  [4]. 

2.  The  working  fluid  and  the  noncondensible  gas  were  treated 
as  thermodynamically  independent  substances. 

3.  Steady  state  conditions  were  assumed. 

4.  The  working  fluid  was  considered  to  be  at  saturation 
conditions  at  all  times. 

5.  Thermal  resistance  of  the  wick-wall  structure,  in  the 
radial  direction,  was  considered  to  be  negligible  compared 
to  the  thermal  resistance  outside  of  the  pipe. 

6.  The  noncondensible  gas  was  treated  as  a  perfect  gas. 
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B.   CONSERVATION  OF  MASS 

1.   Species  'a'  (working  fluid) 

Referring  to  the  differential  volume  element  in 
Figure  2  and  considering  the  system  on  a  unit  time  basis, 

mass  of  'a'  in   =  mass  of  'a'  out 


or. 


p  V  Ac   =   [p  V  +-r-ip   V  )dx]Ac+m  Pvdx  (II.  1) 

a  a         a  a   qx  a  a         a 


which  reduces  to 


0   =  m  Pv  +  Ac  |-(p  V  )  (II. 2) 

a        dx  "^a  a 


Fick's  first  law  of  binary  diffusion,  as  developed 

in  reference  6,  was  used  to  relate  the  mass  flux  of  working 

fluid  (p  V  )  to  the  mass  flux  of  the  mixture  (pv) .   For 
a  a 

the  one-dimensional  case  the  scalar  form  of  Fick's  law  can 
be  written, 

dco 
Pa^a  -  ^a^P^)   =   "P^ab  ^  ^^^'^^ 

where  the  mass  fraction  of  species  'a'  is  defined  as. 


CO   =   p  /p  (II. 4) 

a      a 


and  the  density  of  the  mixture  is  given  by 
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P   =   P^  +  Pb  (II. 5) 


since  species  'a'  and  species  'b'  occupy  the  same  volume  as 
independent  substances. 

Then  equation  II. 2  becomes, 


*a  KS  ^   S'-a^^'  -  lz<P°ab  ^'   =   °  <"-^' 


which  represents  the  conservation  of  mass  for  species  'a'. 
2 .   Species  'b '  (noncondensible  gas) 

A  similar  equation  was  written  for  species  'b'  noting 
that  since  species  'b'  is  a  noncondensible  gas,  m,  =  0 . 
Then, 

^(w,  pv)  -  i-(pD,   -5-2.)   =   0  (II. 7) 

dx  b      dx   ba  dx 

where  the  mass  fraction  of  species  'b'  is  defined  as. 


Wj^   =   Pj^/P  (II. 8) 


3.   Mass  Continuity  of  the  Mixture 

Adding  the  equations  for  species  'a'  and  species  'b' 
(III. 6  and  III. 7)  gave  the  continuity  equation  for  the 
mixture. 


p     ,  ,      doj    d(jo. 


where  D  =  D  ,  =  D,   is  used  for  simplicity  of  notation, 
ab     ba  rr  ^ 
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Recalling  equations  II. 4,  II. 5,  and  II. 8, 


w   +  00,    =1 

a    b 


which  reduced  equation  I I. 9  for  mass  continuity  of  the  mix- 
ture to. 


m  ?^  +  i-(Pv)   =   0  (11.10) 

a  Ac    dx 


C.   CONSERVATION  OF  MOMENTUM 

1.   Mixture  of  Species  'a'  and  Species  'b' 

Conservation  of  linear  momentum  for  the  mixture,  with 
steady  state  conditions,  was  written  for  the  vapor  space 
differential  element  of  Figure  3  as. 


[f   =     /  v(pv  •  dA) 
cs 


where,  dA  is  the  outward  directed  surface  area  element  and, 

->• 

F  is  the  sum  of  all  external  forces  acting  on  the  fluid  in 

the  control  volume.   For  the  one-dimensional  case  F  =  F  . 

X 

The  sum  of  the  body  forces  and  surface  forces  was  seen  from 
Figure  3  to  be, 


If        =      p  g  Ac  dx  -  T  Fv  dx  -  dp  Ac      (11.11) 


and  the  surface  integral,  representing  the  net  efflux  of 
linear  momentum  in  the  axial  direction  from  the  control 
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Figure  3.   Vapor  Space  Differential  Element 
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volume,  was  expanded  to  give: 


/  =   -  pv  Ac  +  ( V  +  dv)  ( pv  +  d ( pv)  ) Ac  +  v  m  Pvdx 

a  a 
cs ,  X 


=  Acpvdv  +  Acvd(pv)  +  v  m  Pvdx  (11.12) 

a  a 


equating  equations  11.11  and  11.12,  and  rearranging  yielded, 


Pv   dp        dv  ,    d(pv)  ,    .   Pv       ,T.^  ,-,. 

pg  -  T  :r -^        =    pv  ^r—  +  V    j,       +  V  HI   tt—  (11.13) 

^     Ac   dx        ax      dx      a  a  Ac 


For  steady  state  conditions  the  noncondensible  gas  remains 
stationary  with  respect  to  the  coordinates  fixed  to  the  heat 
pipe.   Therefore  the  net  flux  of  species  'b'  is  identically 
zero  (i.e.,  Pu^v,  =  0  since  v,  =  0)  .   Noting  that 
pv  =  P^v^  +  p^vj^,  and  Vj^.  E  0, 


pv   =   p  V 
a  a 


or 


_   _p_        V 

a     p        CO 
a        a 


Replacing  v   in  equation  11.13  and  rearranging  resulted  in 

3. 

the  momentum  equation  for  the  mixture, 

m 

dp         d  ,   2,   Pv    'a    Pv  ,^^    ,.v 

d$  =   P5-a3j(PV  )  -^  V  —  -  ^T  (11.14) 

a 
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where  the  viscous  shear  force  t  is  some  function  of  velocity. 

Su 
For  example,  x  =  —  v,  if  Poiseuille  flow  is  assumed. 

2 .   Species  'b' 

Since  the  velocity  of  species  'b'  is  zero,  the  con- 
servation of  momentum  equation  for  species  'b'  could  be 
written  directly  by  inspection  as 


dp^ 


=  Pk  g  (11.15) 


dx      ^b 

A  separate  momentum  equation  for  species  'a',  similar 
in  form  to  equation  11.14  could  be  written.   However  the  addi- 
tional equation  would  not  be  useful  since  it  would  not  be 
linearly  independent  of  the  two  momentum  equations  already 
presented. 

D.   ENERGY  BALANCE  FOR  PIPE  WALL 

In  developing  a  heat  pipe  energy  balance  the  wick  and 
pipe  wall  were  considered  together  as  a  composite.   Recall 
that  thermal  resistance  of  the  wick-wall  structure,  in  the 
radial  direction,  was  assumed  to  be  negligible.   Figure  4 
then  represents  an  appropriate  differential  element  upon 
which  to  develop  the  energy  balance. 

By  the  conservation  of  energy  (referring  to  Figure  4), 


in      out 


then. 
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wick 


T(x) 


m  h^  Pv  dx 
a  fg 


Qc  +  dQc 


pipe  wall 


h  PS  (T  -  T  )  dx 

oo 


X  +  dx 


Figure  4.   Differential  Element  for  Energy  Balance 
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Qc+m^h^  Pvdx  =   Qc  +  dQc  +  h  Ps  (T  -  T_^)  dx  (11.16) 


but,  by  Fourier's  law  of  conduction. 


Qc   =   -  kAp  II  (11.17) 


Incorporating  equation  11.17  into  equation  11.16  and 
rearranging  yields  the  energy  balance  equation. 


d^T  _  hP^(T-T  )  ^  ^^^m    =   0  (II. li 

^^2    k  Ap      oo'     k  Ap    a 


where, 

(a)  the  coefficient  of  thermal  conductivity  (k) 
must  be  selected  to  represent  the  composite 
wick-pipe  wall  structure,  and 

(b)  the  heat  transfer  coefficient,  h,  is  some 
function  of  temperature.   The  correlation 
selected  for  h(T)  depends  upon  the  external 
cooling  conditions. 

E.   EQUATIONS  OF  STATE 
1.   Species  'a' 

Since  saturation  conditions  were  assumed  to  exist 
in  the  working  fluid  at  all  times,  tabulated  saturation 
properties  of  the  working  fluid  can  be  used  to  generate 
functions  for  partial  pressure,  partial  density  and  heat  of 
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vaporization  of  the  working  fluid  as  functions  of  local 
temperature.   Then, 


p    =   Fa(T)  (11.19) 

a 


p    =   Sa(T)  (11.20) 


The  functions  Fa(T) ,  Sa(T) ,  as  well  as  h^  (T)  repre- 
sent  actual  tabulated  saturation  properties  of  the  working 
fluid  which  can  be  readily  used  in  a  computerized  interpo- 
lation scheme.   Alternatively  Fa(T),  Sa(T),  and  h^  (T)  could 
be  approximated  by  polynomials  or  other  anlaytic  functions 
selected  to  fit  the  working  fluid  saturation  curves. 
2 .   Species  'b' 

As  a  perfect  gas, 


Pk   =   Pu  ^^-  T  (11.21) 


'b      ^b 


^ 


3 .   Mixture  Properties 

For  independent  substances. 


P   =   Pa  -^  Pb  (11.22) 


and. 


P   =   P,  +  Pk  (11.23) 

a    b 
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Although  equation  11.22  is  valid  based  upon  Dalton's  law 
of  partial  pressures,  it  is  included  here  for  completeness 
only  since  it  is  not  linearly  independent  of  equations 
11.14  and  11.15  developed  from  mom.entum  considerations. 

F.   SUMMARY  OF  EQUATIONS 

The  following  set  of  ten  independent  equations  in  ten 
unknowns  describes  steady  state  operation  of  a  vertically 
oriented  gas  loaded  heat  pipe. 


*all*a;<p^'    =  °  'i^-"' 


3I  =  pg  -  33^(pv  )  -  _  V  ^  -  ^  T  (11.14) 

a 


^Pb 


=   p,  g  (11.15) 


dx      "b 

d^T   _       PS      ^      ^_^         ^   ^   PV   ^         ^      Q  g 

,  2    k  Ap        °o      k  Ap   a 

b 


00^   =   P  /P  (II. 4) 

a      a 


P   =   P,  +  Pk  (11.23) 

a    b 
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p    =   Fa(T)  (11.19) 

3. 


p    =   Sa(T)  (11.20) 

3. 


The  ten  dependent  variables  (m  ,  v,  p,  p  ,  p,  ,  p,  p  ,  p,  , 

a         a   o      3.        ID 

00  ,  T)  are  all  functions  of  axial  position, 
a 
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III.   TRANSFORMATION  TO  COMPUTER  MODEL 

A.   ELIMINATION  OF  VARIABLES 

The  ultimate  goal  was  to  obtain  two  nondimensionalized 
differential  equations  with  nondimensional  temperature  and 
mass  fraction  of  species  'a'  as  the  sole  dependent  variables. 
By  reducing  the  number  of  equations  and  dependent  variables 
analytically  the  number  of  computations  required  for  each 
numerical  integration  step  is  significantly  reduced.   This 
approach  was  taken  to  enhance  the  numerical  accuracy  and 
stability  of  computer  results  for  the  model.   With  this 
objective  in  mind,  the  set  of  ten  equations  in  ten  unknowns 
summarized  in  Section  II. F  was  reduced  to  two  equations  in 
two  unknowns  by  the  following  steps. 

1.   Total  pressure  (p)  appeared  in  equation  11.14  only  and 
therefore  could  not  be  eliminated  as  a  dependent  variable  by 
back  substitution.   Since  total  pressure  as  a  function  of 
axial  position  was  not  of  particular  interest,  p  was  eliminated 
by  assuming  that  the  axial  pressure  gradient  of  total  pressure 
was  negligible  (dp/dx  =0).   Since  the  axial  pressure  gradient 
of  total  pressure  is  primarily  created  by  the  hydrostatic 
head  of  the  vapor-gas  mixture,  this  common  assumption  is 
reasonable.   This  assumption  means  that  total  pressure  is 
assumed  to  remain  essentially  constant  throughout  the  vapor- 
gas  space  of  the  pipe  for  a  given  operating  power,  gas  load, 
operating  temperature  and  sink  temperature. 
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It  is  significant  to  note  that  although  the  axial  gradi- 
ent of  total  pressure  is  negligible,  axial  variations  in 
the  partial  pressures  of  both  working  fluid  vapor  and  noncon- 
densible  gas  are  significant.   For  example,  using  methanol 
as  the  working  fluid  at  typical  operating  conditions  results 
in  an  axial  partial  pressure  variation  of  several  atmospheres. 

This  assumption  eliminated  equation  11.14  and  the  variable 

P- 

2.   The  following  simple  substitutions  were  then  sequen- 
tially made  in  the  remaining  equations. 

a.   Equation  11.20  was  used  to  eliminate  p  . 

a 


P^   =   Sa(T) 
a 


b.   Equation  11.19  was  used  to  eliminate  p 

a 


p    =   Fa(T) 
a 


c.   Equation  11.21  was  used  to  eliminate  p. 


Pb  =  "b  i^  ' 


d.   Equation  11.23  was  used  to  eliminate  p. 


Pb   =   P  -  ^a 


e.   Equation  II. 4  was  used  to  eliminate  p 


=  ^  =   Sa(T) 

GO  0) 

a       a 
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3.   The  surviving  four  equations  in  four  unknowns  (m  / 

3i 

v,T,ij)    )  were: 

Si 

^aAF  ■"  d3^^^^^^   -  d^^-i;r- -d^^    =    0  (iii.i) 


\E-|^^^   =   0  (III. 2) 

a 

a      D  a 

d^T   PS  h(^_^  ^  ^  ^^^2^m   =   0  (III. 4) 


,2    k  Ap     °o      k  Ap    a 


4.   Equation  III.  2  was  then  used  to  eliminate  m   in 


equations  III.I  and  III.  4. 


Ac  d  ,Sa  V  X 

"^a   =   ■  p^  dSi^-ir-^ 

a 


Equation   III.I   then  became: 


d    ,Sa   Vv  d    ,^  V         d    ,Sa   Da,              « 

-  4-  +  3—  Sa   V     -   ^{^^ ^        =      0 

dx      0)  dx  dx      0)          dx 

a  a 


which  was  integrated  directly  to  yield. 


c    ft          1  \  Sa  D   a     ^ 

S  a  V  (1  -  — )  -  —5 —  =   C 

CO  CO    dx 

a  a 


where  C  is  a  constant  of  integration.   However,  at  x  =  0 
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V   =   0 


and     Sa/co 


then 


do) 


C   =   -  D  p 


o  o  dx 


(III. 5) 


where  x  =  0  was  assumed  to  be  at  the  condenser-reservoir  end 

dco__ 

=  0  is  a  boundary  condition  as 
o 


of  the  pipe.   However, 


dx 


will  be  discussed  later. 


Therefore   C   =   0 


The  remaining  system  of  three  equations  was  then. 


Sa  v(l  -  ^)  -  Sa__D   a   ^   q 

CO       CO    dx 
a      a 


(III. 6) 


^[Sa(—  -  1)#-  T] 
dx     CO       M, 


Sa(—  -  l)g 
a 


(III. 3) 


d  T   Ps  h 


dx 


2    k  Ap 


(T  -  T  ) 


h^   Ac  ,   -, 
fg    d  ,Sa  V 


k  Ap    dx   CO 


(^:^^^^)   =   0 


(III. 7) 


5.   Equation  III.  6  was  finally  solved  for  velocity  and  used 
to  eliminate  v  in  equation  III. 7  to  give  a  system  of  two 
equations  in  T  and  co  . 


|-[Sa(—  -  1)-^  T]   =   Sa(;i-  -  1)  g 

dx     CO       M,  CO       -^ 

a      b  a 


(III. 3) 


j2„    r^   1,  h_   Ac  ,     T^  <^     doo 

d  T  _  Ps  h.  _  m  ^  _   fg     d  r      D  Sa      a 


dx 


2    k  Ap 


k  Ap   dx  CO  (co  -1)   dx 
^      a   a 


=   0 


(III. 


31 


To  solve  these  two  equations  simultaneously  it  was  advan- 
tageous to  expand  and  differentiate  equation  III. 3  and  to 
expand  and  rearrange  equation  III. 8  thereby  generating  a 

system  of  two  second  order  differential  equations  in  'jj 

a 

and  T. 

After  differentiating  equation  III. 3  a  considerable  amount 
of  algebra  led  to  the  following  form: 

dx  a  dx  CO 

a 

^^      oj    ^  ^^     dT    ^^      ^y  ^dx^       ^    ^       2    dx     dx^  ^^     dT    ^^^^ 
a  dT  CO 

a 

00  a  CO 

a  a 

Equation  III. 8  was  manipulated  into  a  similar  convenient 
format  as: 

,2„         d^co      Ac   h^      D  Sa 


dx^  dx^'^P   ^    (-a-l)-a  ^^P  ^ 


Ac  h  doj  dco  ,  , 

+   1^ V    f    "^    ^, ;t^^-D   Sa(-^)  (-ip  +  ^) 

Ap   k    (co   -l)co       dx  dx       co   -1        co 

^  a  a  a  a 

■^    (2^    dT    ■"    ^   "dl^^     d^^  (III. 10) 


Equations  III. 9  and  III. 10  represent  the  final  form  of  the 
desired  equations  in  a  dimensional  format. 
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B.   BOUNDARY  CONDITIONS 

To  numerically  integrate  these  second  order  differential 
equations  four  boundary  conditions  must  be  established  at 
X  =  0.   To  avoid  the  obvious  singularity  when  co   =  1.0  (at 
the  evaporator  end  of  the  condenser)  integration  must  proceed 
from  the  condenser-reservoir  end  toward  the  evaporator. 
Therefore  x  =  0  is  chosen  to  be  the  condenser-reservoir  end 
of  the  heat  pipe. 

Appropriate  boundary  conditions  are  then  as  follows. 


1.      T      =      T  at      X  =   0  (III. 11) 

00 


2.^=0  at      X  =   0  (III. 12) 

dx 

dco 
3.      -^     =0         at      X   =    0  (III. 13) 


4.      0)        =0)        >0        atx=0  (III. 14) 

^  ^o 


The   selection  of  an   initial   value  of  oj      is    somewhat   arbi- 

a 

trary.   However,  the  initial  value  of  co   selected,  in  conjunc- 

a 

tion  with  the  values  of  T   and  T   used  in  the  model,  indirectly 
fix  the  amount  of  noncondensible  gas  and  the  operating  power 
level  of  the  heat  pipe  being  simulated  for  a  particular 
working  fluid  -  noncondensible  gas  pair. 

C.   NONDIMENSIONALIZATION  OF  EQUATIONS 

To  nondimensionalize  equations  III. 9  and  III. 10  the 
following  dimensionless  variables  were  defined. 
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XS   =   x/L  nondimensional  length 

(independent  variable) 

T  -  T 

oo 

Yl   =   = — _  rp        nondimensional  temperature 
e   °o 


where: 


T    =   T(x) 


T   =   prescribed  evaporator  temperature 
T   =   ambient  temperature 


Y2   =   In  03         logarithm  of  the  mass  fraction  of 
^        species  'a' 


For  ease  of  notation  and  computer  programming  the  first 
and  second  derivatives  of  the  two  dependent  variables  were 
represented  as  follows. 


^   d(Yl)    ^ 
^■^      dXS 


=   d(Y2)   ^ 
^^      dXS        ^^^ 


DY3   =   ^^^^^^  =  d^(Yl) 

^^^  dXS^ 

DY4   -   ^(P^2)  -  d^(Y2) 

"    ^^^  "  dXS^ 


Each  function  of  temperature  was  nondimensionalized  with 
respect  to  the  function  value  at  evaporator  temperature  (T  ) 
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The  following  nondimensionalized  functions  of  temperature 
were  then  developed  for  direct  substitution  into  equations 
III. 9  and  III. 10 . 


Gl   = 


D 


e 


nondimensional  binary  mass 
diffusion  coefficient 


G2   =   ^^  ^''e-J^l 


dT    D(T  ) 
e 


nondimensional  temperature 
gradient  of  the  binary  mass 
diffusion  coefficient 


G3   = 


12. 


fg   e 


nondimensional  heat  of 
vaporization 


G4   = 


G5   = 


Sa 


Sa(T  ) 
e 


d(Sa)  ^'^e"'^=«^ 


dT     Sa(T  ) 
e 


r,6  -  d'(Sa)  "e-^J 


nondimensional  working 
fluid  density 


nondimensional  temperature 
gradient  of  working  fluid 
density 


nondimensional  curvature  of 
working  fluid  saturation 
density  curve 


G7   = 


h(T  ) 
e 


nondimensional  convection 
heat  transfer  coefficient 


Using  the  nondimensional  variables  and  functions  of 
temperature  defined  above,  equations  III. 9  and  III. 10  were 
transformed  into  the  following  nondimensionalized  system  of 
two  equations  in  two  dimensionless  dependent  variables  (Yl 
and  Y2) ,  and  their  first  and  second  derivatives  with  respect 
to  axial  (dimensionless)  length. 
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DY3{  (1  -  Exp(Y2)  )  (  (K1Y1+1)G5  +  K1G4)  }  -  DY4 (  (KlYl+l)G4) 


(l-Exp(Y2) ) {Y3(Y3(2K1G5+(K1Y1+1)G6)+K2G5) } 


+   Y4{2Y3( (K1Y1+1)G5+K1G4)    -  G4 ( (KlYl+1) Y4-K2)  }  (III. 15) 


DY3(Exp(Y2)-l)  -  DY4(K3G1G3G4)   =   (Exp ( Y2) -1) K4G7Y1 
+  G1G3Y4(K3(G4(^||1  -  ^Expa2Hl)^  +G5Y3))  (III. 16) 

where  Kl  through  K4  are  dimensionless  constants: 


(III. 17) 


T  -  T 
Kl   =  {-^ ) 

00 


00    -^  /=• 


Ac  h^  (Tp)  D(T  )  Sa(T  ) 
i^r>   =  (  fg   ^ e e_ 

^   k  Ap  (T  -  T  ) 
^   e   =° 

Ps  h(T  )  L^ 

K4  =   ( — r—TT^ ) 

k  Ap 


In  terms  of  the  nondimensional  variables,  the  boundary 
conditions  (equations  III. 11  through  III. 14)  became: 

1.     Yl   =   0  at   XS   =   0  (III. 18) 


2.     Y2   =   Y2   <  0        at   XS   =   0  (III. 19) 

o 
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3.     Y3   =   DYl   =   0  at   XS   =   0         (III. 20) 


4.     Y4   =   DY2   =   0  at   XS   =   0         (III. 21) 


where  Y2  =  In  co  .   As  noted  previously  with  the  dimensional 

form  of  this  boundary  condition  (equation  III. 14),  the  initial 

value  of  CO   selected  is  arbitrary.   However  the  initial  value 
a  ^ 

of  CO  ,  in  conjunction  with  the  specified  values  of  T   and  T  , 
a.  e      °° 

indirectly  specifies  the  quantity  of  noncondensible  gas  and 
operating  power  level  of  the  heat  pipe.   Since  0  <  co   <  1.0, 

3. 

the  corresponding  value  of  Y2  must  be  negative. 
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IV.   DISCUSSION  AND  RECOMMENDATIONS 

A.   DISCUSSION 

The  Continuous  System  Modeling  Program  (CSMP) ,  developed 
by  the  IBM  Company  [5],  is  a  straightforward  computer  package 
available  to  numerically  solve  equations  III. 15  and  III. 16. 
The  CSMP  package  includes  a  number  of  optional  numerical 
integration  methods  including  a  variable  step  fourth-order 
Runge-Kutta  method. 

Although  equations  III. 15  and  III. 16  are  highly  nonlinear 
in  the  dependent  variables  Yl  and  Y2,  this  system  of  equations 
can  be  considered  linear  in  the  second  derivatives  of  Yl  and 
Y2 .   Therefore  these  equations  can  be  written  in  the  following 
form: 

(COEFl) (DY3)  +  (C0EF2) (DY4)   =   Fl  (IV. 1) 

(C0EF3) (DY3)  +  (C0EF4) (DY4)   =   F2  (IV. 2) 

where, 

COEFl   =   { (l-Exp(Y2) ) ( (K1Y1+1)G5+K1G4) }  (IV. 3) 

C0EF2   =   -( (K1Y1+1)G4)  (IV. 4) 
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Fl      =       (l-Exp(Y2) ) {Y3(Y3(2K1G5+(K1Y1+1)G6)+K2G5) } 

+    Y4{2Y3(  (K1Y1+1)G5+K1G4)-G4( (KlYl+l)Y4  -  K2) } 

(IV. 5) 

C0EF3      =       (Exp(Y2)-l)  (IV. 6) 

C0EF4      =      -(K3G1G3G4)  (IV. 7) 

F2      =       (Exp(Y2)-l)K4  G7  Yl  +G1  G3  Y4(K3(G4 


. fG2  Y3     _    Y4Exp(Y2)  rcy3^ ^  . ^V    8^ 

^^r^         (Exp(Y2)-l)    ^   G5Y3))  (IV. 8) 


Using  these  definitions,  the  nondimensionalized  equations 
representing  steady  state  heat  pipe  operation  may  be  numerically 
integrated  using  the  variable  step  fourth-order  Runge-Kutta 
option  of  the  CSMP  routine  as  follows. 


Yl      =         /       (DYl)dXS      =         /      ^4|^  dXS  (IV.  9) 

0  0 


Y2      =         /       (DY2)dXS      =         /      ^i^l^    dXS  (IV. 10) 

0  0 


1  12 

Y3      =         /       (DY3)dXS      =         /      ^    ^^^^    dXS  (IV.  11) 

0  0         dXS 
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1  12 

Y4   =    /   (DY4)dXS   =    /   ^  ^^^^  dXS  (IV. 12) 


0  0     dXS^ 


where. 


DYl   =   Y3  (IV. 13) 


DY2   =   Y4  (IV. 14) 


and  by  Cramer's  rule, 


DY3   =    (F1C0EF4  -  F2  C0EF2) 

(COEFl  C0EF4  -  C0EF2  C0EF3  )  vxv.x^j; 


nv4   =    (F2  COEFl  -  Fl  COEF  3) 

(COEFl  C0EF4  -  C0EF2  C0EF3  )  Vxv.xd; 


This  model  can  be  used  to  demonstrate  gravitational  effects 
on  the  vapor-gas  interface  region  by  providing  axial  profiles 
of  temperature  and  the  mass  concentration  of  the  working 
fluid.   When  the  molecular  weight  of  the  noncondensible  gas 
is  less  than  that  of  the  working  fluid,  the  axial  profiles 
will  shew  a  fairly  small  interface  region  between  vapor  and 
gas.   However,  when  the  molecular  weight  of  the  noncondensible 
gas  is  substantially  greater  than  that  of  the  working  fluid, 
the  profiles  will  show  a  distorted  (spread  out)  interface 
region. 

Qualitative  results  can  be  obtained  from  this  model  by 
using  additional  simplifying  assumptions  and  built-in  CSMP 
functions  as  follows: 

40 


1.   Nominal  constant  values  can  be  selected  to  represent 
the  convection  heat  transfer  coefficient  and  the  binary  mass 
diffusion  coefficient.   For  example, 

2 

h(T)   =   1.0  BTU/HR/FT    for  free  convection 

D(T)   =   0.4  FT  /HR   for   Krypton- Methanol 


or 


2 

D(T)       =      1.9    FT   /HR      for      Helium-Methanol 


2.   To  generate  values  of  working  fluid  density  (Sa)  and 
heat  of  vaporization  (h_  )  at  any  temperature,  a  CSMP  function 
generator  can  be  used  to  provide  quadratic  interpolation  of 
tabulated  saturation  properties.   Tabulated  thermodynamic 
properties  of  methanol  are  available  in  reference  7.   The 

CSMP  function  "SLOPE"  can  be  used  to  generate  values  of 

2       2 
d(Sa)/dT  and  d  (Sa)/dT  as  functions  of  temperature. 

The  temperature  profile  generated  by  this  model  for  a 

given  vertical  heat  pipe  using  a  particular  working  fluid- 

noncondensible  gas  pair,  with  given  input  values  of  T  , 

T^ ,  and  oo  (0)  ,  represents  heat  pipe  operation  at  a  certain 

power  level  and  total  pressure  with  a  specific  load  of  non- 

condensible  gas.   The  power  level,  total  pressure,  and  mass 

of  noncondensible  gas  represented  could  be  computed  as 

follows. 
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Mass   of   Gas    ( lb    )       =     ^  D   ^   L      /      — (l-oj    )  dXS  (IV. 17) 

m  4      V  ^-^       CO  a 

0  a 


1 
Power   Level    (BTU/HR)       =      Pv  L       /      m      h^      dXS  (IV. 18) 

0     a   tg 


Sa(T  )  R  T 

00  ' 


Total  Pressure  (psia)   =   Fa(T  )  + t^tt — (l-oo  (0)) 


oj  (0)  '    a   '  M, 
a  D 

(IV. 19) 

B .   RECOMMENDATIONS 

Based  upon  preliminary  work  using  CSMP,  numerical  results 
from  this  model  could  be  improved  by  the  following  measures. 

1.  Analytical  expressions  are  recommended  in  lieu  of  the 
special  CSMP  functions  to  represent  working  fluid  saturation 
properties  as  functions  of  temperature.  Expressions  for  the 
first  and  second  derivatives  of  density  as  functions  of  tem- 
perature could  then  be  found  analytically  and  used  directly. 
The  resulting  expressions  should  yield  more  accurate  values, 
particularly  for  the  second  derivative  of  density,  than  are 
achievable  through  multiple  use  of  the  CSMP  "SLOPE"  function 
with  tabulated  saturation  values. 

2.  Since  the  binary  mass  diffusion  coefficient  is  a  fairly 
strong  function  of  temperature,  quantitative  results  in  the 
interface  region  could  be  improved  considerably  by  use  of 

an  empirical  expression  for  D(T)  instead  of  a  nominal  constant 
value.  Reference  8  contains  appropriate  empirical  expressions 
for  binary  mass  diffusion  coefficients. 


42 


3.  Quantitative  results  could  also  be  improved  somewhat 
by  using  an  appropriate  empirical  expression  for  the  convec- 
tion heat  transfer  coefficient  as  a  function  of  temperature 
instead  of  a  nominal  constant  value.   Reference  9  provides 

a  variety  of  expressions  for  h(T)  depending  upon  the  convec- 
tive  condition  being  simulated. 

4.  With  the  improvements  noted  above,  this  model  could  be 
used  in  an  iterative  loop  to  obtain  preliminary  quantitative 
design  information  concerning  vertically  oriented  terrestrial 
heat  pipes  when  the  molecular  weights  of  the  operating  sub- 
stances are  substantially  different. 

For  example,  if  the  working  fluid- noncondensible  gas  pair 
is  specified,  along  with  the  required  heat  source  (evaporator) 
temperature  range,  ambient  conditions,  and  the  maximum  safe 
total  pressure  for  structural  integrity  of  the  pipe,  then 
equation  IV.  19  could  be  used  to  determine  the  value  of  oj  (0) 
to  use  for  boundary  condition  III. 19.   Equation  IV. 17  could 
then  be  used  to  numerically  compute  the  required  mass  of  non- 
condensible gas  to  load  depending  upon  the  heat  pipe  dimen- 
sions.  Equation  IV. 18  would  then  yield  the  maximum  heat 
rejection  rate  obtainable  for  these  conditions.   If  this  heat 
rejection  rate  is  not  adequate  for  the  design  requirements, 
the  designer  could  vary  any  of  a  number  of  parameters,  such 
as  operating  temperature,  pipe  dimensions,  total  pressure 
allowed,  substances  used,  etc.,  until  an  adequate  design  is 
found  by  iteration. 


43 


LIST  OF  REFERENCES 


NASA  Contract  Report  2018,  Theory  and  Design  of  Variable 
Conductance  Heat  Pipes,  by  B.  D.  Marcus,  April  19  72. 

Bienert,  W.,  Heat  Pipes  for  Temperature  Control,  Proc. 
4th  Intersociety  Energy  Conversion  Engineering  Conference, 
Washington,  D.  C,  1969,  p.  1033. 

Edwards,  D.  K.,  and  Marcus,  B.  D.,  "Heat  and  Mass  Transfer 
in  the  Vicinity  of  Vapor-Gas  Fronts  in  Gas-Loaded  Heat 
Pipe",  Journal  of  Heat  Transfer,  Trans.  ASME  Series  C, 
Vol.  94,  No.  2,    May  1972,  pp.  155-162. 

Kelleher,  M.  D.,  and  Batts ,  W.  H.,  Effects  of  Gravity  on 
Gas-Loaded  Variable  Conductance  Heat  Pipes,  Proc.  2nd 
International  Heat  Pipe  Conference,  Bologna,  Italy,  1976, 
p.  225  . 

Continuous  System  Modeling  Program  III  (CSMP  III) , 
Program  Reference  Manual  SH  19-7001-3,  Program  Number 
5734-XS9,  4th  ed. ,  IBM  Corporation,  Data  Processing 
Division,  White  Plains,  New  York,  December  1975. 

Bird,  R.  B.,  Stewart,  W.  E.,  and  Lightfoot,  E.  N., 
Transport  Phenomena,  1st  ed,  p.  502,  John  Wiley  and  Sons, 
Inc.,  New  York,  1960. 

Smith,  J.  M. ,  "Thermodynamic  Properties  of  Methyl  Alcohol", 
Chemical  Engineering  Progress,  v.  44,  pp.  521-528,  July, 
1948. 

Reid,  R.  C,  Prausnitz,  J.  M.  ,  Sherwood,  T.  K.,  The 
Properties  of  Gases  and  Liquids,  3rd  ed.,  pp.  544-601, 
McGraw-Hill  Book  Company,  New  York,  19  77. 

Holman,  J.  P.  Heat  Transfer,  4th  ed.,  McGraw-Hill  Book 
Company,  New  York,  19  76. 


44 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  California   939  40 

3.  Department  Chairman,  Code  69Mx  1 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

4.  Associate  Professor  M.  Kelleher,  Code  69Kk     1 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

5.  CDR  David  C.  Kane,  USN  1 
4466  Jefferson  St. 

Napa,  California   94558 


45 


Kane  ^9147!     I 

Analytical  treat 
^!f  Of  gravit^ 
^^fects  on  ^a/i      . 
variflhi^  loaded 

heat  t       ^°"^^^*ance 
i^eat  pipes. 


Thesis  191  471 

K135T  Kane 

c.l         Analytical  treat- 
ment of  gravity 
effects  on  gas  loaded 
variable  conductance 
heat  pipes. 


I  treatment  of  gravity  effects 


3  2768  002  11396  1 

DUDLEY  KNOX  LIBRARY 


Vvi)!  vX 


':r:--^:\'M 


■    V"Cs'LVv>-AJ 


-iBiii^lliiti 


mm 


